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Abstract 

A catastrophe model for the onset of fast magnetic reconnection is presented that suggests 
why plasma systems with magnetic free energy remain apparently stable for long times and then 
suddenly release their energy. For a given set of plasma parameters there are generally two stable 
reconnection solutions: a slow (Sweet-Parker) solution and a fast (Alfvenic) Hall reconnection 
solution. Below a critical resistivity the slow solution disappears and fast reconnection dominates. 
Scaling arguments predicting the two solutions and the critical resistivity are confirmed with two- 
fluid simulations. 

PACS numbers: 52.35.Vd, 52.65.-y 
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Explosive events in plasmas, such as solar eruptions and the sawtooth crash in labora- 
tory fusion devices, are driven by magnetic reconnection. Understanding the mechanism 
facilitating fast reconnection in high temperature plasma systems has been a long-standing 
challenge. Sweet-Parker (SP) reconnection Q] is far too slow to explain observations and 
Petschek reconnection requires the invocation of anomalous resistivity, a phenomenon that 
is at best only poorly understood. A new paradigm has emerged in recent years in which 



dispersive whistler and kinetic_A 
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fven waves facilitate fast reconnection by setting up the 

n □ 

, |5( . Magnetospheric satellite observations [6fl and recent 



laboratory experiments support this new paradigm [3] . 

It is not sufficient, however, to merely explain how fast reconnection can occur. If re- 
connection were always fast, magnetic stress could never build up in plasma systems such 
as the solar corona and the explosive release of magnetic energy seen in nature and the 
laboratory would never occur. It is critical, therefore, to explain why fast reconnection does 
not always take place. We show that there are generally two reconnection solutions for a 
given set of parameters: slow reconnection as predicted by Sweet and Parker; and fast col- 
lisionless reconnection facilitated by coupling to dispersive waves in the dissipation region 
(Hall reconnection). Below a critical resistivity the slow solution disappears. The emerging 
picture, therefore, is that slow reconnection can dominate the dynamics of a system for long 
periods of time but the resulting rate of reconnection is so slow that external forces can 
continue build up magnetic stresses. When the resistivity rj drops below a critical value (or 
the available free energy crosses a threshold) the system abruptly transitions to fast recon- 
nection and is manifest as a magnetic explosion. Such a model complements earlier ideas 



that the onset of solar flares, for example, results from the loss of MHD equilibrium 
or more complex "breakout" models 10| | . 



A rather simple argument can be made to motivate why magnetic reconnection is bistable, 
i.e., has two solutions for a given set of parameters. The SP solution is valid provided the 
half width of the current layer 5 exceeds the relevant kinetic scale lengths 
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where L is the half length of the SP current sheet, di = c/u P i is the ion inertial length, p s 
is the ion Larmor radius, u P i is the ion plasma frequency, and ca is the Alfven speed. The 
Alfven speed is to be evaluated immediately upstream of the current layer. Therefore, if the 



system is undergoing SP reconnection and the resistivity is lowered, SP reconnection will 
continue as long as Eq. (JTJ) is satisfied. 

Conversely, fast reconnection is validprovided the dispersive (whistler or kinetic Alfven) 
waves that drive kinetic reconnection |5[ are not dissipated. We restrict our discussion to 
whistler waves, generated by the Hall term. The dispersion relation for resistive whistler 
waves is u = k 2 CA<ii — ik 2 rjc 2 /Air. Since both terms scale like k 2 , dissipation can only be 
neglected if it is small enough at all spatial scales, that is 
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X" < CAd,. (2) 

This can also be written as v e i <C Q C e, where u e i = r]ne 2 /m e is the electron-ion collision fre- 
quency and Q ce = eB/m e c is the electron cyclotron frequency, a condition which is typically 
easily satisfied in nature. Therefore, if the system is undergoing Hall reconnection and the 
resistivity is increased, it will stay in the Hall configuration as long as Eq. (J2J) is satisfied. 
If the resistivity is an intermediate value such that both Eqs. (JIJ and (JSJ) are satisfied, then 
either solution is accessible and the system is bistable. 

We now present estimates of the slow-to-fast (rj s f) and fast-to-slow (r// s ) resistive transi- 
tion boundaries of the bistable regime. We estimate T] s f by setting the left and right hand 
sides of Eq. (JTJ equal using d{ as the relevant kinetic scale length for Hall physics: 
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To estimate r/f s , we perform a Sweet-Parker type scaling analysis [2j that is more precise 
than the argument used to motivate bistability in Eq. (J2J). Resistive effects are negligible if 
the outward magnetic diffusion across the electron current sheet, rjc 2 / '4n5 2 , is less than the 
inward convection, Vi n /5, where Vi n is the flow speed into the electron current layer. For Hall 
reconnection, numerical simulations JjJ have shown that 5 scales like the electron inertial 
length d e = c/oj pe , where uj pe = \J 'Anne 2 j 'm e is the electron plasma frequency, and the inflow 
speed scales like t> in ~ O.lcAe, where CAe is the electron Alfven speed based on the magnetic 
field immediately upstream of the electron current layer. The critical resistivity r]f s is found 
by equating the two: 



2 

C ^in n CAe 



or, using c Ae d e = c A d h 
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~ O.lcxdi, (4) 
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FIG. 1: Normalized reconnection rate, E' , as a function of island width, u>, for the two sets of 
simulations described in the text. The vertical dotted lines show when the added effects were 
enabled. Note that the final parameters of the two solid line simulations are identical. 

where ca is evaluated upstream of the electron current layer. This is consistent with Eq. ©, 
but more precise since the geometry of the layer is included. Note that r]f s is independent 
of system size and electron mass and is enormous for most physical systems. Equation (j3J) 
suggests that once Hall reconnection onsets, resistive effects are unlikely to influence the 
dynamics. The ratio of Eqs. (jSJ) and (@J) gives f] s f/Vfs ~ 10dj/L <C 1, which is small because 
di <C L for most systems of physical interest. Thus, bistabilty is present over an enormous 
range of resistivity. 

The predictions of this model are amenable to tests using numerical simulations. We 

n 

use the two-fluid code, f3d, a massively parallel code described elsewhere |12l| . to perform 
two-dimensional simulations in a slab geometry of size L x x L y . The initial equilibrium is 
two Harris sheets, B = itB tanh[(y ± L y /4)/w ] with wq = 2di, in a double tearing mode 
configuration with periodic boundary conditions in all directions. The ions are initially 
stationary and initial pressure balance is enforced by a non-uniform density. For simplicity, 
we treat an isothermal plasma. A coherent perturbation to induce reconnection is seeded 
over the equilibrium magnetic field. The resistivity i] is constant and uniform. We use small 
fourth-order dissipation, oc 774V 4 with r/4 = 2 x 10 -5 , in all of the equations to damp noise 
at the grid scale. 

The computational domain must be chosen large enough to have a discernible separation 
of scales between the SP and Hall reconnection rates, but with high enough resolution to 
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distinguish the electron inertial scale. We find that a computational domain of L x x L y = 
409. 6dj x 204.8g?j, with a resolution of As = Ay = O.ldi and an electron to ion mass 
ratio of m e = rrii/25 (i.e., d e = 0.2<iA is sufficient. Since the rate of Hall reconnection is 
insensitive to the electron mass , we do not expect our results to depend on our 

particular choice of m e . For this computational domain, we can estimate r] s f and rjf g . In 
evaluating Eq. (jSJ), we use L ~ L x /A = 102.4c?j. Normalizing lengths to di and velocities to 



cao = Bq/ \J Ann^mi, where no is the initial density far from the sheet, we obtain 

To evaluate Eq. (0J), we use the value of B ~ 0.3Bq upstream of the electron current layer 
measured in the simulations to evaluate ca, so 

c 2 
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A larger system would produce a greater separation between rj' s f and r]'j s and would be closer 
to the parameters of real systems but would be more computationally challenging. 

To demonstrate bistability of reconnection with a resistivity in the intermediate region 
v'sf < v' < Vf s i we perform two related sets of simulations. First, we show that a system 
undergoing Hall reconnection with a resistivity below rj'j s continues Hall reconnection for 
any value of resistivity below this value. We start with a benchmark collisionless (r/ = 0) 
Hall-MHD simulation that is run from t — until the rate of reconnection is steady. The 
normalized reconnection rate E' = cE/BqCao is shown as a function of island width w as the 
thick solid line in Fig. ^ The reconnection rate is calculated as the time rate of change of 
magnetic flux between the X-line and O-line. The rate of reconnection jumps to E' ~ 0.06 
by the time the island width is 10di, after which it is remains steady. When w ~ 35di, we 
enable a resistivity of rj' = 0.015 (which lies between the predicted values of r]' s j- and 77^) 
and continue the simulation until most of the available magnetic flux has been reconnected. 
For comparison, the thick dashed line shows the reconnection rate when we maintain rf = 0. 
Clearly, the reconnection rate remains nearly unchanged after the inclusion of the resistivity. 



For the second set of simulations, we want to show that a system undergoing SP recon- 
nection continues to reconnect at the lower rate for any value of resistivity exceeding r] s f. 
Our computational approach is to disable the Hall and electron inertia terms and evolve the 
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FIG. 2: (Color online) Out of plane current density, J z , for late times from the two solid lines 
of Fig. ^ The top plot corresponds to the thick solid line (Hall reconnection) . The bottom plot 
corresponds to the thin solid line (SP reconnection). 

resistive system with a resistivity that exceeds rj s f. We then re-enable the Hall and electron 
inertia terms and continue to advance the full equations. This benchmark simulation is 
performed with rj' = 0.015 (the same value of resistivity as in the run shown in the thick 
solid line in Fig. [1}, and the reconnection rate is again plotted in Fig. Q as the thin solid 
line. The reconnection rate remains stationary with E' ~ 0.01, a factor of six slower than 
the Hall case even with the Hall and electron inertia terms enabled. For comparison, the 
thin dashed line in Fig. ^ shows the reconnection rate for a system in which the Hall term 
is not enabled. Thus, the Hall and the electron-inertial terms do not impact the rate of SP 
reconnection for these parameters. 

The out of plane current density, J z , is shown at late time in Fig. |2] for the runs cor- 
responding to the two solid curves in Fig. ^ The top plot corresponds to the thick solid 
curve. The current sheet is short and opens wide, as is expected in Hall reconnection 
16l Il7t Il8l |19| The bottom plot corresponds to the thin solid curve. The current 
sheet is long and thin as is expected from the SP theory of resistive reconnection [2^, 21, 22^ . 



Since the same equations govern the two sets of data and the value of the resistivity is the 
same, we conclude that the system is bistable. 

To complete the mapping of the two reconnection solutions, we vary the resistivities of the 
benchmark Hall and SP reconnection solutions of Fig. ^ For the case of Hall reconnection, 
corresponding to the thick solid line in Fig. ^ we change rf from 0.0 to 0.010, 0.013, 0.015, 



6 



, 0.08 

a) 

0.06 

bj 0.04 

0.02 
0.00 

0.000 0.005 0.010 0.015 0.020 0.025 0.030 
V' 

b) 1 A ^~ 

1-2 r . 
1.0 - • 

A 0.8 r 
*= 0.6 r 
0.4 r 

0.2 : Q • « • o • o o o ° 

0.0 : 

0.000 0.005 0.010 0.015 0.020 0.025 0.030 

7] 

FIG. 3: (a) Steady state normalized reconnection rate, E', as a function of normalized resistivity, 
7/ for runs analogous to those in Fig. [I\&s described in the text, (b) Current sheet width, 5, as a 
function of 7/ for the simulations in (a). 

0.0175, 0.020, 0.0225, 0.025 and 0.030 when w ~ 35^. For the case of SP reconnection, 
corresponding to the thin solid line in Fig. ^ we change 7/ from 0.015 to 0.003, 0.007, 
0.009, 0.011, 0.013, 0.0175, 0.020, 0.0225, 0.025 and 0.030 when w ~ 50rf, (after the Hall 
and electron-inertial terms have been re-enabled). The asymptotic reconnection rate is 
computed as the time averaged reconnection rate once transients have died away. 

The results are plotted in Fig. Eta), with the states starting from Hall reconnection plotted 
as open circles and the states starting from SP plotted as closed circles. The closed circles 
reveal that the disappearance of the SP solution occurs abruptly, with T]' s f between 0.011 
and 0.013. The open circles reveal the disappearance of the Hall reconnection configuration, 
with t]'j s between 0.020 and 0.0225. The error bars are due to random fluctuations in the 
reconnection rate. The plot is reminiscent of what one would expect of a bifurcation diagram 
for a system with a cusp catastrophe. 

Thus, the numerical simulations confirm that magnetic reconnection is bistable over a 
range of resistivity consistent with the scaling law predictions of rj' s f ~ 0.01 and rj'j- s ~ 0.03. 
The asymptotic steady state current sheet width 5, calculated as the half width at half 
maximum of J z (y) at the X-line, is plotted in Fig. [2Kb) for each of the runs. As predicted 
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FIG. 4: Normalized reconnection rate, E', as a function of island width, w, for the simulation 
which is started at rf = 0.015, reduced to 0.007, then increased back to 0.015. 

by Eq. (JIJ), the steady state SP current sheet width 5 is of order when the resistive 
reconnection solution ceases to exist, as is shown by the closed circles of Fig. Efb). 

We emphasize that the results presented in Fig. El though generated by a specific nu- 
merical procedure, are not sensitive to the details of this procedure. To demonstrate this, 
we show that the key feature of Fig. EH the boundary where the slow reconnection solution 
disappears can be reproduced through a hysteresis-like procedure: in the simulation corre- 
sponding to the thin solid line in Fig. Q we first lower the resistivity from rf = 0.015 to 
rj' = 0.007 (when the island width is about w ~ 50dj). As expected from Fig. Eland shown 
in Fig. HJ the transition from SP to Hall reconnection occurs. We then raise the resistivity 
back to rj' = 0.015 (the original value) when the island width is about w ~ 68dj. As can be 
seen in Fig. El fast reconnection continues, showing that the system can be in either of two 
steady states for the same set of parameters. 

Having verified that reconnection is bistable with Hall MHD simulations, we return to the 
onset problem. We can compare the critical resistivity (or temperature, using the classical 
Spitzer formula) with observations of onset in physical systems such as solar eruptions and 
sawtooth crashes. For solar flares, n ~ 10 10 cm -3 , L ~ 10 9 cm and B ~ 100 G j23], so from 
Eq. (13 , rj s f ~ 10~ 16 s in cgs units, corresponding to a temperature of 10 2 eV ~ 10 6 K. This 
is in excellent agreement with the coronal temperature. 

For the sawtooth crash, the relevant kinetic scale is p s instead of dj because of the presence 
of a large guide field. If we assume that Eq. (JTJ) can be carried over to the guide field case, 
we can make a comparison for sawtooth crash onset. Typical parameters for sawteeth in the 
DIII-D tokamak |24| are B v ~ 2 T, T e ~ 2.0 keV,r s ~ 20 cm,n ~ 10 14 cm" 3 and zT cff ~ 2, 
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and L ~ r s 8, where 9 ~ 60 ° is the angular extent of the current layer 25|. For bean-shaped 
flux surfaces, the helical field strength in the plasma core is B ~ 100 G, so using Eq. Q with 
p s in place of dj yields r] s f ~ 10~ 16 s, corresponding to T ~ 2 x 10 2 eV. This temperature is 
an order of magnitude too small. However, the inclusion of diamagnetic effects |26j], which 
are known to slow reconnection, should improve agreement. In future work we will explore 
whether Eq. (0) does hold in the presence of a guide field, for which p s is the relevant kinetic 
length scale. 

The effect of collisionality on the reconnection rate was recently explored in the Mag- 
netic Reconnection Experiment (MRX) j^. A sharp increase in the reconnection rate was 
observed at low collisionality, consistent with the qualitative picture presented here. Data 
for the current sheet width is unavailable, so quantitative comparisons are not possible at 
this time. Further, in this experiment fast reconnection has been correlated with magnetic 
turbulence localized in the reconnection layer j^. Since the present simulations are limited 
to 2-D we can not address the development of this turbulence and how it might impact our 
conclusions. We surmise that our conclusions will not be strongly changed as long as the 
turbulence does not broaden the reconnection layer beyond the scale length di. 
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